Nishimori point in the 2D ±J random-bond Ising model 
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We study the universality class of the Nishimori point in the 2D ±J random-bond Ising model 
by means of the numerical transfer-matrix method. Using the domain-wall free-energy, we locate 
the position of the fixed point along the Nishimori line at the critical concentration value p c = 
0.1094 ± 0.0002 and estimate v = 1.33 ± 0.03. Then, we obtain the exponents for the moments of 
the spin-spin correlation functions as well as the value for the central charge c = 0.464 ± 0.004. The 
main qualitative result is the fact that percolation is now excluded as a candidate for describing the 
universality class of this fixed point. 
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In the past years, the subject of disordered systems has 
known a huge renewal of interest in the condensed mat- 
ter and statistical mechanics community. Among these 
disordered models, two-dimensional systems are of a par- 
ticular interest. Since the discovery of the unitary series 
of conformal field theory (CFT) in 1984 Q, exact values 
for the exponents of many well known models of statis- 
tical mechanics have been given. However, an equivalent 
classification for universality classes of such systems in 
the presence of impurities is still missing. A first big step 
towards a more general classification has been done re- 
cently as a random matrices classification 0. The data 
for critical exponents in most of the experimental rele- 
vant fixed points for impure system is however still not 
available. 

The Ising model on a square lattice is one of the most 
popular two-dimensional systems. It is specified by the 
energy of a spin configuration 
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where the sum is over all bonds and the coupling con- 
stants Ji.j are bond dependent. We consider here the 
Jij = ±1 Random-Bond Ising Model (RBIM) with the 
following probability distribution: 

P(J iyj ) = P 8{J^ - 1) + (1 - p)S{J i , j + 1) . (2) 

Note that with these conventions, the pure model (p = 0) 
is characterized by Jij — — 1 and thus has a ferromag- 
netic groundstate. 

The RBIM is similar to other relevant disordered mod- 
els such as the Chalker-Coddington random network 
model which was proposed originally in the context of 
the quantum Hall effect plateau transition [pi. However, 
it is important to stress that these systems have a differ- 
ent phase diagram and therefore their fixed points have 
no reason to be in the same universality class H. 

The topology of the phase diagram of the RBIM de- 
pends crucially on the type of disorder one considers. An 



instructive example is provided by a disorder having only 
two possible values for the bonds with equal signs and 
probabilities. It is by now well established || that the 
only non-trivial fixed points are located at the extrema of 
the boundary of the ferromagnetic phase, corresponding 
to the pure Ising fixed point and a zero temperature fixed 
point which turns out to be in the percolation universal- 
ity class. It is interesting to notice that percolation is 
also the universality class of the so-called spin quantum 
Hall model Q , another random network model. 
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FIG. 1. Phase diagram of the two-dimensional ±J ran- 
dom-bond Ising model. 

When the distribution contains also bonds with differ- 
ent signs (like in (§)), the situation is more subtle. For 
a certain class of probability distributions, Nishimori has 
shown that a so-called 'Nishimori' line exists where many 
properties can be calculated exactly @. For the proba- 
bility distribution (0), this line is given by 



1-p 
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with (3 = l/T. On the Nishimori line, the internal en- 
ergy can be calculated exactly and an upper bound can 
be given for the specific heat. Also of interest is an equal- 
ity of the moments of the spin correlation functions (see 
below) . Nishimori has further proven inequalities for the 
correlation functions which yield important constraints 
on the topology of the phase diagram which is shown in 
Fig. for the ±J RBIM [§. Since the Nishimori line is 
also invariant under Renormalization Group (RG) trans- 
formations [[t0| , the intersection of the Nishimori line and 
the Ferro-Para transition line must be a fixed point. This 
so-called Nishimori point (TV) corresponds to a new uni- 
versality class belonging precisely to the family of strong 
disorder fixed points. The bold line in the phase dia- 
gram Fig. [| is the phase boundary between the ferromag- 
netic and paramagnetic regions. At zero temperature the 
model has a spin glass phase |11|. The three non-trivial 
fixed points along the bold line are the pure Ising fixed 
point, the Nishimori point at the crossing with the Nishi- 
mori line (dotted line) and the zero temperature point, 
separating the ferro and the spin glass phases. The prop- 
erties of the latter point are still mostly unknown. This 
point will be the subject of future investigations [Ti| . 

In the last years, many numerical and analytical efforts 
have been made in order to identify the universality class 
of the Nishimori point. Very recently, an analytic ap- 
proach suggested that it is governed by an Osp(2n+l\2n) 
symmetry (and maybe Osp(2n + 2\2n) but unfortu- 
nately, the classification of CFTs with such symmetries 
is still missing. From a numerical point of view, there is 
by now a long list of results |13-17|. An important ob- 



servation is that all the numerical results for the critical 
exponents tend to suggest that this point is in the per- 
colation universality class. Because of the importance of 
the statistical model on its own and its relevance for un- 
derstanding the plateau transition in the quantum Hall 
effect, it is crucial to elucidate the similarity to percola- 
tion and the relation to the super-symmetric CFT pro- 
posed in the literature. 

In this letter we provide results of extensive numerical 
transfer-matrix calculations of the Nishimori point with 
the binary distribution (^) for bonds on the square lat- 
tice. We use the domain-wall free-energy to accurately 
locate the critical concentration of disorder p c and to es- 
timate the exponent v. Then, we analyze the spin corre- 
lation functions and the scaling of the free energy, giving 
accurate and novel results for the magnetic exponent r\ 
and central charge c. Apart from improving the identi- 
fication of this universality class (providing in particu- 
lar values for the central charge which have never been 
measured before), our main result is that percolation is 
excluded as a possible candidate for describing this fixed 
point. 

First, we use the free energy of a domain wall 
to locate the critical point. For a strip of width L the 
domain- wall free-energy dh is defined as |l8| 
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where fjf* is the free energy per site of a strip of width 

L with periodic boundary conditions and the corre- 
sponding one with antiperiodic boundary conditions, di, 
is an observable which can be used directly to study the 
RG flow under scale transformations. In particular, it is 
constant at a fixed point. 



We have computed f£ 
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ploying a standard transfer matrix technique with sparse 
matrix factorization (see, e.g., Jn|) on strips of length 
N = 10 6 . Since randomness is strong, care must be taken 
to reduce fluctuations even if the free energies are self- 
averaging. Therefore, we have fixed the concentration 



(p) 



and 



of bonds p globally on a sample and computed f L 
f[ a ' on the same sample. Still, one needs around 1000 
to 4000 samples of L x 10 6 strips to obtain sufficiently 
small error bars for L < 12. Even on modern computers 
this needs an amount of CPU time which precludes the 
analysis of wider strips. However, since we are looking 
for a fixed point, no crossover effects are expected and it 
is legitimate to use small system sizes. 
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FIG. 2. Domain-wall free-energy. The inset shows the raw 
data and the main panel the scaling collapse with p c = 0.1094 
and v — 1.33. The symbols are for L — 8 (filled boxes), L — 9 
(open boxes), L = 10 (*), L = 11 (x) and L = 12 (+). Error 
bars are much smaller than the size of symbols. 

The inset of Fig. || shows dz,(p) along the Nishimori 
line (J3J) in the vicinity of the critical concentration p c . A 
finite-size estimate for p c is given by the crossing points 
^ii(Pc) = dL 2 {Pc)- After extrapolation to an infinitely 
wide strip (details will be given elsewhere p2|), one ob- 
tains 



p c = 0.1094 ±0.0002 



(5) 



This estimate improves upon the accuracy of earlier esti- 
mates |p^-|l7|. It agrees perfectly with the transfer ma- 
trix computations fl4|,[l7| while we find a slightly smaller 



2 



value of p c than [|5|,|6| . We would like to mention that 
(||) is confirmed by standard Monte Carlo simulations on 
systems up to 32 x 32 sites Q - the present estimate is 
just more accurate. 

One can extract also the correlation length exponent 
v from d,L if one assumes the scaling form 

d L (p- Pc )=d((p- Pc )L 1 /^ . (6) 

Again, we omit details |l2| and quote just the final result 
v= 1.33 ±0.03. (7) 



demonstrates that follows 
) with the parameters (0) and 



The main panel of Fig. 
indeed the scaling form 

®' 

The result (|7|) is in complete agreement with v = 1.32± 
0.08 obtained by high-temperature series p"5j] a s well as 
the value v = 4/3 for percolation (see, e.g., ]20[). 

Another important quantity is the magnetic exponent 
77. This exponent can be measured, for example, by com- 
puting spin-spin correlation functions. As we mentioned 
in the introduction, all along the Nishimori line the mo- 
ments of these correlation functions are equal two by two: 

[(S{x liyi )S{x 2 ,y 2 )) 2k - 1 ] = [(S{x uVl )S{x 2 ,y 2 )) 2k ] (8) 

for any integer k. Here [■ • ■] stands for the average over 
the disorder. Assume now that the correlation functions 
(||) decay algebraically on a plane and define by x, y the 
coordinates on the infinite cylinder of circumference L, 
with x G [1,L] and y g] — oo,+oo[. Using a conformal 
mapping, one infers then the following behavior of the 
correlation functions on the cylinder: 



i(S(x u y)S(x 2 ,y)) n ]K 
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(9) 



For a pure system, one has rj n = n x 77. On the other 
hand, in the case of percolation over Ising clusters, it is 
easy to see that the moments of spin correlation functions 
are all equal (and not only two by two). Then, if the 
Nishimori point is in the percolation universality class, 
the exponents for the correlation functions in (||) should 
collapse to a unique value rj n = 77 at the critical point. 

In order to verify this we have calculated the spin- 
spin correlation functions on cylinders of width L and 
length 400 x L, {i.e. with the length >• L) for L up to 
20. We have checked that for width L = 12, lattice and 
finite length corrections are of order 1%. One example 
of these correlation functions can be seen in Fig. || (with 
x\ = y = and x 2 = x) on a doubly logarithmic scale. 
One observes that the correlation functions nicely obey 
the power law (^|), thus verifying both the correct location 
of the critical point as well as the functional form of the 
spin-spin correlation function in a finite strip. 
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FIG. 3. Moments of the spin-spin correlation function for 
p = 0.1095 and L = 20. We only show the odd moments: 
n = 1 (+), n = 3(x),n = 5(*) and n = 7 (open boxes). 
Error bars are smaller than the size of the symbols. The 
values of the exponents are given in (|lo|). 

We can then fit the exponent by studying the depen- 
dence with distance of the correlation functions (Q) or 
by studying the dependence with L for the fixed location 
x = L/2. The first method has proven to give smaller 
error bars and we obtain for the family of exponents r\ n 
for p c = 0.1095 and L = 20: 
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= 0.3354, 



(10) 



with relative errors at most of the order of 1%. 
One immediately notices two things: 

i) The value for r/i differs considerably from the value 
of percolation 77 = 5/24 w 0.2083 (see, e.g., |o)), 

ii) the exponents for higher moments are also consider- 
ably different from 771 which is also clear from inspection 
of Fig. |. 

These results are compatible with the behavior of the 
magnetic susceptibility that will be presented elsewhere 
jl2| . We have also calculated estimates for the exponents 
assuming different values for p Cl namely, p c = 0.109 and 
0.110 and the results are still distinct from the ones of 
percolation (we obtain 771 = 0.180(1) for p = 0.109 and 
i]i = 0.190(1) for p — 0.110). Moreover, it is only in 
the region very close to p — 0.1095 that we obtain a 
stable estimate for r/i as we increase the width L of the 
lattices. One can then conclude from the exponents con- 
trolling the algebraic decay of the correlation functions 
that the Nishimori point is not in the percolation univer- 
sality class. 

Finally, we discuss the (effective) central charge c. It 
characterizes the number of gapless degrees of freedom at 
the critical point and appears as the universal coefficient 



3 



of the first finite-size correction to the free energy for 
periodic boundary conditions 21 
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(11) 



The leading term f<£ is not universal and does indeed 
change already when we modify the conventions for the 
model (Q). There are higher-order finite-size corrections 
to the free energy including terms of the form L~ 4 . The 
central charge is an important quantity identifying the 
CFT description JlJ of a fixed point. One has c — 1/2 for 
the critical point of the pure Ising model, but it has not 
been determined yet for the Nishimori point. 

In the process of computing dr, we have also obtained 
estimates of f L p ^ for different values of p. One can either 
fit these values for exactly by ( |ll| ) ignoring further 
corrections in which case the data for the smallest values 
of L should not be used. Or one includes a correction 
term of the form L -4 which improves the convergence 
with system size. These two approaches yield consistent 
estimates for a given p. In addition, one can test that 
the result does not change significantly if other higher- 
order corrections are added. It should also be noted that 
the sensitivity of the estimates for c with respect to the 
location of p c is negligible in comparison with the errors 
coming from the finite-size analysis. The final result is 
that the following is a safe estimate for c at the Nishimori 
point of the ± J RBIM (more details will be given in [fl2"l|): 



0.464 ±0.004. 



(12) 



Assuming again the universality class of percolation, we 
would expect the value for percolation in the Ising model 
c = 5v fj" 2 w 0.4777 [§. Even if our result (|l|) is close 
to this value, it can still be distinguished safely from 
percolation. This finding is one more argument that the 
Nishimori point is not in the universality class of percola- 
tion, at least the one expected from Ising clusters. Notice 
also that the central charge does not rely on a choice of 
observables, and this argument can then be considered 
as the most general one. 

The results presented in this letter provide new insight 
into the 2D ±J RBIM and the Nishimori point. The 
main result is that, according to the magnetic exponent 
and the behavior of the higher moments of the corre- 
lation functions, the universality class is different from 
the one of percolation, at least considering the Ising spin 
variables as fundamental observables. This conclusion is 
supported by the central charge which is independent of 
the choice of observables and was measured here for the 
first time. These results and a detailed study of the zero 
temperature fixed point will provide a general and com- 
plete control of the RBIM, which one could consider as 
the simplest, but most fundamental model for disordered 
systems in 2D. 
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